Entropy 2010, 12, 1799-1832; doi:10.3390/el2071799 



OPEN ACCESS 



entropy 

ISSN 1099-4300 

www.mdpi.com/journal/entropy 

Article 

Roofs and Convexity 

Arinin Uhlmaim 

Institute for Theoretical Physics, University of Leipzig, PB 100 920, D-04009 Leipzig, Germany; 
E-Mail: armin.uhlmann@t-online.de 

Received: 5 June 2010; in revised form: 5 July 2010 / Accepted: 9 July 2010 / 
Published: 20 July 2010 



Abstract: Convex roof extensions are widely used to create entanglement measures in 
quantum information theory. The aim of the article is to present some tools which could 
be helpful for their treatment. Sections 2 and 3 introduce into the subject. It follows 
descriptions of the Wootters' method, of the "subtraction procedure", and examples on how 
to use symmetries. 
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1. Introduction 

One often is in the position to know a quantity for pure states of a quantum system, say g, without 
a definite meaning in classical physics. Then one looks for a method to extend g to all mixed states, 
see [1-4]. Let G denote such a possible extension. 

A reasonable approach is certainly to extend g "as linearly as possible" or, more correctly, "as affine 
as possible". The mathematical model for such a demand is a "roof" as presented in Section 2. 

As it is difficult, to say the least, to imagine higher dimensional geometry, let us look at an elementary 
example, the real qubits. They fill a disk bounded by a circle. The circle represents the pure states. 
Considering g{'K), with vr pure, as the height of a wall at the point vr. An extension G of g io the disk 
provides a covering of the ground floor. Let G{ijj) be a point of that covering. To be a roof in the sense of 
Section 2, we like to have: There is either a straight line or a plane coincident with the roof point G{ijj) 
and resting on the wall at two or at three points tt^, g{pij). If the line or plane is parallel to the ground 
floor, it is called flat. An uncomfortable feature is the rather large arbitrariness: There are plenty of roof 
extensions allowing, however, as a bonus, much room for playing. Under rather weak assumptions, there 
is a minimal as well as a maximal roof extension (Proposition 2.1). 
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Physically stronger motivated seem extensions G of g which are either convex or concave. 
Heuristically, the convex ones try to suppress the "classical noise" or else the "classical information" 
by the tendency to attach lower values to the states "far" from the pure ones. The concave extension 
stress the "classical part", possessing lower values in the vicinity of the pure states. 

It is a well known fact that there are maximal convex and minimal concave extensions. If 
it exists, the minimal roof extension is equal to the maximal convex extension and this facilitate 
calculations sometimes. 

The maximal convex extension of g can be gained by 



with pure states tcj, [5]. In quantum information theory this is nowadays a common procedure to 
define entanglement measures. Its first appearance is in the important paper [6] by C. H. Bennett, 
D. DiVincenzo, J. Smolin, and W. Wootters, and has been called entanglement of formation and is 
denoted by Eiuj). These authors considered states w of a bipartite quantum system. The role of g plays 
S'(Trfe7r), the entropy S of the partial trace Tr^, as a function of a pure bipartite states vr. 

In the present paper a trace preserving positive map T from the states of a quantum system into another 
one is called a stochastic map. A channel is a completely positive stochastic map. If in the definition of 
entanglement of formation, E, a stochastic channel T is used in place of the partial trace Tr;,, we call the 
resulting quantity the entanglement ofT and denote it by E^- For a channel, Et is equivalent (in many 
ways) to the restriction of E onto a face of a bipartite state space of sufficiently high dimension, [7]. In 
this way, T is seen as a sub-channel of the partial trace: The tensor product is partitioned into subspaces 
on which the sub-channels are defined. A 1-qubit channel, for example, can be represented by a 2 x m 
bipartite quantum systems as the restriction of the partial trace over the larger dimensional part onto the 
density operators supported by a suitable 2-dimensional subspace. 

Clearly, Et is of interest in its own as part of the x*-function x* = St — Et with {St) (w) = S{T{uj)) 
the maximum of which is the Holevo capacity [8] of T. 

Since the importance of E has been realized in [6], several other measures of similar structure 
have been introduced and discussed, replacing the von Neumann entropy S by another function on the 
state space of the output system. The perhaps mostly discussed examples are the "tangle" tt and the 
"concurrence" Ct- The connection between them is Cr(7r)^ = rr(7r) for all pure states vr. Sometimes 
one also needs a minimal concave extension ("entanglement of assistance"). All these quantities have 
been defined by global variational problems of type (1), (2). 

A further, even earlier source for the said procedure roots in the problem of defining a "quantum 
dynamic entropy", generalizing the Kolmogorov-Sinai one. In the approach [9] of A. Connes, 
H. Narnhofer, and W. Thirring several similar global variational problems wait to be solved. One of 
them is the search for the convex roof defined by g{uj) = S{Dn{co)), S again von Neumann's entropy 
and Dn the diagonal map, setting all off-diagonal entries of a matrix to zero. The problem initiated the 
paper [10] of F. Benatti, H. Narnhofer, and A. Uhlmann and further ones. 




(1) 



where the inf is running through all convex decompositions 




(2) 
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It is true that complete solutions for these extensions are only known in the lowest non-trivial 
dimensions. In the present paper general tools are presented to facilitate the treatment of convex or 
concave roofs. The general aspects are mainly in the Sections 2 and 3. It includes the foliation of 
the input state space into "convex leaves" onto which the roof becomes affine and, in particular nice 
situations, even constant, [10,11]. 

Section 4 is devoted to the Wootters' way [12-14] of presenting the entanglement of formation 
explicitly, see also [15]. 

Section 5 shows a more recent way to compute the concurrence of a 1-qubit stochastic map by 
a substraction procedure. For the concurrence of 2 x m quantum systems it allows to compute the 
concurrence of any rank two quantum state. An elegant way to do so was opened by Hildebrand, [16,17], 
using the so-called "S-lemma"of Yakobovich, see [18]. Another one has been chosen by Hellmund and 
Uhlmann, [19], who use the description of general 1-qubit maps given by Gorini and Sudershan, [20]. 

For the tangles the pioneering work goes back to Coffman, Kundu, and Wootters, [21], who already 
remarked that optimal decompositions of length two should be sufficient in the 1-qubit case. That the 
"substraction procedure" works can be read off from a paper of Osborne and Verstraete, [22]. Here 
we describe analytical results for axial symmetric channels. There are also results for the 3-tangle 
roof problem, [25]. 

Most results, if not numerically, are found by the help of symmetries. Besides the already quoted 
ones, an essential step has been done by K. G. VoUbrecht and R. F. Werner, [26], and B. M .Terhal and 
K. G. VoUbrecht [27]. Meanwhile it became a very large domain of research, exceeding the frame of 
the present paper. Hence, in Section 6, only some aspects, connected with maximal symmetric states, 
are touched. 

Some notations: We use "H for Hilbert spaces, B{l-L) for its algebra of operators, fi('H) for the set 
of density operators supported by "H. We also say "state" for "density operator" and "state space" for 
Vl{1-L). As a convex set, Vl{l-L) is embedded in Herm('H), the real linear space of Hermitian operators. 
The symbol Vl is also used for a general compact convex sets in a real linear space. Following [3], the 
extremal points of a convex set will be called "pure" ones mostly. They are usually symbolized by the 
letter tt. We follow [28] in using ri{x) = — xlogx and S{uj) = Ti r]{uj), the von Neumann entropy. 

2. Roofs, Roof Extensions 

We are now going to give an exact meaning to the word "roof". For this purpose we assume G to be 
a real valued function on a compact convex set Q, contained in a finite dimensional real linear space. 

Definition 2.1a: Roof points 

w G is called a roof point ofG, if there is at least one extremal convex decomposition 





(3) 



such that 




(4) 
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If this takes place, we call the decomposition (3) optimal with respect to G or, equivalently, G-optimal. 
The number of terms in (3) with pj 7^ is the length of the decomposition. 

Definition 2.1b: Flat roof points 

A roof point of G is called y?a?, if there is a G-optimal decomposition (3) fulfilling 

G{uj) = G(7ri) = G(7r2) = . . . (5) 

i.e. all the values G{nj) are equal one to another [29]. 

Let /(x) be a real function defined on the range of G. The main merit of a flat roof point a; of G is 
the simple fact, that it remains a flat roof point for /(G): f{G{TTj)) = f{G{uj)) for all j and (4) remains 
true for /(G). In other words, the flat points of p — G(p) are also flat roof points of p — ^ /(G(p)). 

Definition 2.2: Roofs, flat roofs 

A real function G on 1] is a roof if every cu E ^lis a roof point of G. 
G is called aflat roof if all roof points are flat ones. 

In important applications the point of view is a bit different: A function n ^ g(n) is given on the set 
j^pure puj.g states, or, more generally, on the set of extremal states of an arbitrary compact convex set. 
Then one asks for meaningful extensions G which coincides with g on ^1^^^^. In such a generality the 
problem is too arbitrarily posed and one asks for restrictions to such an extension. Remarkable ones are 
the roof extensions. They interpolate the values of g "as linearly (or as affine) as possible". 

Definition 2.3: Roof extensions 

Let 7r — )■ g{TT) be a real function on fiP^^'^'^. A roof G is called a roof extension of g, if G(7r) = g{'n) for 
pure states. 

Now observe the following simple fact: The maximum (the minimum) 

max{Gi, . . . , G„} respectively min{Gi, . . . , G„} 

of finitely many roof extensions of (7 is a roof again. 

Indeed, assume for roof extensions Gi, . . . , G„ of (7 and uj E Vt the value Gi(w) is not less than the 
other values Gj{uj). Then one selects a Gi-optimal decomposition for max{Gi, . . . , G„}. 

The reasoning above fails for infinitely many roof extensions. The proof of the following is postponed 
to that of proposition 3.5. 

Proposition 2.1 

Let ghea real continuous function on QP^^^{1-L) and G a roof extension of g. For all u E nCH) there 
exist optimal decompositions (3), (4) the length of which does not exceed (dim?^)^ + 1. 
There is a minimal and a maximal roof extension of g. 
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2.1. Examples 

Roof extensions exist in abundance. To see this and also the difficulties, we may get into, let us 
consider a few examples showing some typical constructions. We remain within the state space Vt{T-i), 
dim v. = d and we start with d = 2. 

Example 2.1: A Bloch ball construction 

Seen from convex analysis, the space of all 1-qubit density operators is a 3-dimensional ball, the 
Bloch ball. ilP^^"^ is the Bloch sphere, the surface of the Bloch ball. 

Now assume there is a function g on the Bloch sphere and we like to find roof extensions G of g. 
Particular nice ones are gained as following: We take a bundle of straight lines such that every point u 
of the Bloch ball is coincident with exactly one line, say L^, of the bundle. If to is not pure, then L^^ hits 
the Bloch sphere at exactly two points, say tti and tc2. Hence, a; is a convex combination of them. Now 
we define 

G{uj) = pgi-Ki) + (1 - p)g{TT2) if w = piTi + (1 - p)vr2 (6) 

Because there is just one line from our bundle going through to, we get a well-defined roof. 

Let us now specify our example by choosing a bundle of parallel lines. It belongs exactly one main 
axis of the Bloch ball to the bundle. We may assume that it is the xs-axis with respect to the Bloch 
coordinates xi, X2, of a general Pauli representation 

u = ^{1 + xiai + X2(r2 + x^as) (7) 

The line is now fixed by the values xi and X2, letting the third Bloch coordinate arbitrary. 
crosses the Bloch Sphere at the pure states with the Bloch coordinates Xi, X2, y-i = ±(1 — — x'^Y^'^. 
Hence 

TT± = -{l + Xiai+X2(T2±[l- xj-xlY^'^as), u = pn+ + (1 - p)n_ (8) 
with < j9 < 1. Expressing now g in terms of Bloch coordinates, we obtain the roof 



G{uj) = G{xi, X2) = pgixi, X2, +^1 - xj- xl) + (1 - p)g{xi, X2, -\J I- x\- xl) (9) 

The example is further specified by choosing for g the Shannon entropy of the diagonal elements of 
oj. Then G becomes the roof 



Er(.) = H( '^^'l''--''- , '-^'-/^-''- ) (10) 

Tiuj) denotes the diagonal part of uj and H the Shannon entropy of the diagonal elements. That this 
is the solution for the entanglement of the diagonal channel, has been suggested by Levitin, [30] and 
Thirring [31]. (10) provides a roof extension as described above. Indeed, it is a flat roof depending 
onl — x\ — x\, hence a function of the Euclidean distance ^Jx\ + xi of uj from the a;3-axis. It is 

(0|u;|l) = ^1^^, xl + xl=A\{Q,u\l)\^ (11) 
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Therefore, the concurrence of the diagonal map for qubits can be written 

CT{uj) = \{<d\uj\l)\ (12) 

It is easily to be seen that C is a flat convex roof. It possesses quite large domains where it is just 
affine: Cut the Bloch ball with a plane containing the Xa-axis. We get a disk, cut by the Xa-axis into two 
half-disks. The concurrence is affine on every such half-disk. 

Indeed, if T is a 1-qubit channel with two different pure fix points, one can observe a similar 
phenomenon. (The disks will be cut by the axis through the fix points.) 

Example 2.2 

There is a constructions, similar to the previous one: Given any function /(xs) on the xs-axis. We 
construct a roof by 

A{uj) = f{x^), X3 = (0|a;|0)-(l|a;|l) (13) 

More generally, given a bundle of planes so that every point of the Bloch ball is coincident with one 
and only one of these planes. A function, constant on planes, is a roof. This can be achieved by choosing 
a function / on the Xs-axis and attaching to every plane the value of / at its crossing with this axis. 

Similar constructions are possible with higher dimensional balls (ellipsoids). However, on higher 
dimensional state spaces things are essential more complicated due to the subtle structure of the set of 
their pure states. 

Example 2.3: Affine functions on the state space 
An affine function on Herm('H) is of the form 

X ^1{X) := a + TiXA (14) 
with an Hermitian operator A and a real number a. It is sometimes useful to write (14) in the form 

/(X) = TrX5, B = A+'^1 

(Jj 

The function (14) is trivially a roof: If we have any extremal decomposition (3) one immediately gets 
/(p) = YliPjK'^j)- Thus, every extremal decomposition is optimal. 
We easily conclude that with G also G + Z is a roof. 

A bit more tricky is the assertion that every function (14) is diflat roof on the state space Vl{1-L). 
We prove this by induction to the dimension of the Hilbert space. We start with the qubit case. We 
can choose a basis in H such that A = a + a'a^ and 

u = -{1 + xiai + xsCTs) 

It suffices to prove the assertion for A = a^, resulting in /(X) = Tr a^X = X3. With an unimodular 
number e we consider the pure states 
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so that our affine function is constant at the segment pu^ + (1 — p)uj^. The segment contains uj for 



Varying e we see the following: To every affine function / as in (14) there is an axis through the 
Bloch ball such that I is constant on every plane perpendicular to that axis. We can say something more: 
Given u E Vt and two affine functions, Ij, j = 1,2, we use the constants lj{uj) = Cj to define two 
planes by lj{X) = cj. Because the two planes contain to, they intersect along a line, say L^, containing 
u). The intersection of L^^ with the Bloch sphere provides two pure states which define a flat optimal 
decomposition of cu. 

Proposition 2.2 

Affine functions on ^7('H), dim?/ < oo, are flat roofs. Given two affine functions and a density 
operator cu, there is a common flat optimal decomposition of tu. 

Proof. The 2-dimensional case has already been proved. Let lj{X) = aj + Tr AjX, j = 1,2, denote two 
affine functional and choose to E VL{1-Ld)- Assume the assertion is true for dim?/ < d.liu is from the 
boundary of VL(1-Ld) we are done by our induction hypothesis. Otherwise we consider the linear subspace 
C of Herm('Hd) orthogonal to Ai and A2, i.e., of all Y satisfying Ti AjY = 0, j = 1,2. Consider 
the affine space = u + C. It contains u and it is lj{Y) = lj{co) for all Y E C^^. The intersection 
K = Q n Cui is compact and its extremal points belong to the boundary of VL{1-Ld). Therefore, we 
get a decomposition 



and, furthermore, = Ijiojj) because of Uj E K. However, the support of any ujj is of dimension 
less than d. Therefore, there are extremal decompositions 



P = 1/2. 





k 



for all j and all k with pure states vr^fc. Thus 




3 



k 



is a flat optimal decomposition of u for li as well for I2. Thus the proposition is proved. 



□ 



The following is a corollary. 



Proposition 2.3 



Let f{xi,X2) be a function defined on the range of two affine functions, li and 12- Then 



(15) 



is a flat roof. 
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Example 2.4: An application to the diagonal map 

We like to apply the proposition above to the diagonal map 

X ^ Dn{X) = diag(X) (16) 

which is a channel on Vl{1-L). diag(X) denotes the diagonal part of X obtained by replacing all 
off-diagonal elements of X by zeros. Proposition 2.3 proves that the von Neumann entropy 

uj ^ S{D^{u)) = Y,v{{3\uj\j)) (17) 

is a flat roof for n = 2, 3. For = 3 one replaces the third diagonal element, 0:33, by 1 — xn — X22 to see 
that proposition 2.3 suffices to verify the assertion. 

We like to prove the flatness of (17) for all dimensions. This will be done in example 3.2 later on. 

Example 3 

A general way to obtain roof extensions in state spaces is presented next. It allows for further 
modifications, but it is difficult to be controlled explicitly. 

Let (7 be a continuous real function on QP^^'^ and u E Vt. Every basis \'4ji), . . . , \ilJd) gives rise to a 
decomposition 

a; = ^v^|^,-)(V^,-|v^ (18) 

of positive rank one operators. After normalization of the rank one operators we get extremal convex 
decomposition of uj, 

vm^^ ^^^^ 

In the following definitions we vary over all bases. 

Gi{u) = mm > — , , I I , , — (20) 

''''''' ^0 (^il^l^j) 

where pj is as in (19). Because continuity of g is assumed, and the set of bases is compact, the minimum 
in (20) will be attained. Thus, Gi is a roof extension of g. 
One gets another roof extension by 

G2{uj) = max > pjg{ — / / — ) (21) 
Remark: One cannot guaranty the roof property without continuity of g. 

(21) and (20) are global optimization problems for which there are no explicit expressions known in 
most cases. There are, however, algorithms to approximate them numerically. 

3. Roofs and Convexity 

We start with some definitions and elementary, mostly well-known statements around maximal convex 
and minimal concave extensions. The second subsection is concerned with sufficient conditions for the 
existence of convex and concave roofs and their main properties. 
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3.1. Convex and Concave Extensions 

Let be a compact convex set and g a real function defined on the set VP'^^'^ of its extremal points. 

Definition 3.1: Convex (concave) extensions 

A convex function G is called a convex (respectively concave) extension of g if G is convex 
(respectively concave) and coincides on VP™^ with g. 

Between any convex, concave or roof extension of a function g the inequalities 

^convex ^ ^roof ^ ^concave (22) 

are valid. Indeed, with an optimal decomposition (3) for C^""^ we get 

^^"^(u;) = ^p^.G"^°°f(7r,) = J]pjG™"(7r,) 

because the extensions coincide on VP^'^^. The right sum cannot be smaller than G'^°^^^^{ijj) by convexity. 
Similarly one argues in the concave case. 

The proof of (22) is valid pointwise, leading to 

Proposition 3.1 

Let (^convex ^ convex, (^concave ^ concave, and G any extension of (7. If w is a roof point of G, then 

^convex < G{uj) < G'"'"'''™ (w) (23) 

Let G^^o"™^ and (^^oncavc -j^ j-j^^ above proposition and assume in addition equality in (23), 

(^convox(^^^ ^ G™"=^^^(cj). For all convex combination 

(jj = ^^pjWj, CJj G ^, Pj > 0, (*) 

we conclude 

j^p-G^'^iuj) > G™"(u;) = G'=°^^'^™(w) > J]j9jG^°°"^"°(wj) (24) 

because of (22) the conclusion is G^°'^™^(wj) = G'=°°=^^<=(wj) for all j. 

To express this finding we need some standard terminology. A subset K C ^l is called a. face of if 
from CO G K and (*)it necessarily follows coj E K for all j. Faces are convex subsets of Vl. If not only Vt 
but also QP^^^ is compact then faces are compact and any face K is convexly generated by QP"^'' fl K. 
The intersection of faces is either empty or a face again. The smallest face containing uj will be called 
u-face of f2 and denoted by face(^[f2]. If is a face and to not a point of the boundary of K, then K is 
the u-face of il. 

Now we express the finding above by 
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Proposition 3.2 

Coincide a convex and a concave extension of (7 at w G Vt, then they coincide on the w-face of Vt. 

The least upper bound of a set of convex functions is convex again. Hence, given g on VP^^^, there is a 
unique largest convex extension of g from QP^^'^ to fi. Similarly there exists a smallest concave extension 
oi g. It is convenient to introduce an extra notation for these extensions: 

Definition 3.2: g'^ and g'^ 

Let (7 be a real function on r2P"''^. We denote by g^ the largest convex and by g"^ the smallest concave 
extension oi gioVt, 

g^ = largest convex extension of g from n^^^'^ to , 
g'^ = smallest concave extension of g from ^V^'^^ to . 

We also write G = (or G = C^) if G is the largest convex (or the smallest concave) extension of 
the restriction of G onto il^^^'^. 

Proposition 3.3 

Let G be an extension of g and to one of its roof points. If G is convex, then G{uj) = g^{uj)- If G is 
concave, then G{uj) = g'^{uj). 

If a convex (resp. concave) roof extension of g exists, then it is unique. Because G is convex, G < g^ . 
Because w is a roof point, (23) asserts G{uj) > g'^ico). 

Is there a convex extension at all for a given g 1 1f there is one then there is also a largest one, i.e. g^ 
exists. The answer to the question is affirmative and has been given in [5] by a variational characterization 
which is well known in quantum information theory as a recipe to construct entanglement measures: 

g'^iuj) = MJ^pM^j) (25) 
g'^iu) = sup^pjgiTTj) (26) 

where the "inf", respectively "sup", is running over all extremal convex decompositions 

uj = J2Pj'^j^ TTjG^^P"''' (27) 

of u. 

Indeed, if G is an extension of g which is convex, the right side of (25) must be always larger than 
G{uj). On the other hand, given ui and U2, one can find decompositions (27) for them differing an 
arbitrary small amount e > o from g^{coi) respectively g'~'{uj2). They may be composed from the pure 
states TTjj and probabilities Pij, i = 1,2. Then 

2e + p(?^(u;i) + (1 -p)(?^(u;2) > J^m.i^^K,,) + ^(1 - p)(7(7r2,fc) 

and this is not smaller than g'^{puji + (1 — p)uj2)- Because e can be arbitrary near to zero, g'~' is convex. 
The concave case can be settled by a similar reasoning or by 

-9'' = i-9r, -9'' = {-9r (28) 
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Remark: Hulls of functions 

The convex hull of G is the largest convex function which is smaller than G. The concave hull of G 
is the smallest concave function which is larger than G. In (25) and (26) one uses the values of G at the 
pure states only. Because the hull construction must respect values on the whole of Vt, there are more 
constraints to be fulfilled. 

The expressions (25) and (26) are similarly structured as those of the convex and the concave hulls of 
a function G onVt. One mimics the proofs and gets 

conv[G](w) = iYii^pjG{ujj) (29) 
coYi.c[G]{uj) = sup y^^pjG^Uj) (30) 

where, as in (24), one has to run through all convex combinations 

oj = ''^^pjUj, Uj G Q, Pj > 

Obviously, conv[G] < G^ and conc[G] > G'^. This quite simple reasoning provides also 

Proposition 3.4 

If G is a concave or a roof extension of g then = conv[G]. 
If G is a convex or a roof extension of g then ^f*^ = conc[G] . 

As a matter of fact one can do similar hull constructions with any subset of ^7 which convexly 
generates il. This has been emphasized in [26]. 

3.2. Convex and Concave Roofs 

If there is a convex roof extension of g, then it is equal to g'~'. There is a sufficient condition to 
guaranty the roof property of g'~' and of g^. 

Proposition 3.5 

Let be a convex set. Assume both, and QP"''^, are compact and g continuous on f^P"''^. Then g'~' 
and g'^ are roofs. According to proposition 3.3 they are the minimal respectively maximal roof extensions 
of ^. 

Remember that ^l^^^^ is compact if f2 = n{'H) and "H is finite dimensional. The requirement of 
continuity of g is often satisfied in physically motivated applications, though not always. A counter 
example is the Schmidt number in bipartite quantum systems. In this case it is not known whether g^ 
and are roofs. Nevertheless, the assumptions needed are rather weak ones. We met them already in 
proposition 2.1. 

The proof will be "constructive" in a certain sense. It is arranged to sharpen "theorem 1" in [10]: If 
we know an optimal decomposition with pure states tti, 7r2, . . ., then every convex combination of them 
is optimal. In particular, the (convex or concave) roof is affine on the convex set generated by the pure 
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states TTi, 7r2, . . .. Restricted to this set, the graph of G is a piece of an affine space. The whole graph 
of G appears as composed of affine pieces [32]. If one would know a covering of Vt by these "convex 
leaves", one could compute from the values of g at VP^^^. Things are similar for g^. 

To start proving propositions 3.5 and 2.1 let us repeat the assumptions, is a compact convex set in 
a real linear space £ of finite dimension, the set QP"^^ of all pure {i.e., extremal) points of Vt is compact, 
(yf is a real continuous function on QP"""^. The dimension of f2 as a set in C is denoted by n. It is the 
dimension of the affine space generated by Vt. 

Remark: The space Vl{1-L) of density operators is embedded in Herm('H). The latter is of dimension d? if 
dimH = d. The affine space generated hy Vl{'H) is the hyperplane of Hermitian operators of trace one. 
The dimension oi Vl{1-L) h n = — 1. 

We enlarge C to the linear space £' = £ © IR. Its elements, X © A, will be written in vector form 
{X, A} with two components, X E C and A G H. We need the set 

E = {TX,g{Tx)}, TTGr^P"'^^ (31) 

is a compact set by our assumption. Hence its convex hull, denoted by VL[g], is a compact convex set. 
The set of extremal points of VL\g] is E. (If one of the elements of E would be a convex combination of 
the others, the same would be true for the corresponding pure states, contradicting our assumptions.) 

Choose uo eVI and consider in £ © IR the straight line consisting of the points {to, A}, A G IR. The 
line intersects with VL[g] along a compact segment 

eVL[g] ^ \o{uj) <\<\i{u) (32) 

A satisfies (32) if and only if there is an extremal decomposition 

{a;, A} = $^p,{7r„^(7r,)} = {uj, J^pM^^) >' ^ (33) 

Therefore, 

g'^iu) = Xo{uj) < A < Xiioo) = g'^{uj) (34) 

and there exist extremal decompositions of u with equality in (25) respectively (26). Therefore, g^ and 
g'^ are roofs and proposition 3.5 has been proved. 

If G is a roof extension of g, the point {u, G{uj)} is contained in ^l[g]. Hence it can be represented by 
a convex combination of elements from E. As the dimension of ^l[g] is n + 1, there are, by a theorem of 
Caratheodory, pure convex decomposition of length n + 2. This proves proposition 2.1: For f2 = ^1{'H) 
it follows n + 2 = - 1) + 2. 

On the other hand, a point {cu, g'~'{uj)} belongs to a face of the boundary of ^l[g]. Its dimension cannot 
exceed n. Thus there are, again by Caratheodory, pure decompositions of length n + 1. For ^Kji) this 
gives an achievable length d"^, an often used fact [33]. 

Definition 3.2: convex leaves 

Let G be a real function on n. A subset K C ^lis called a convex leaf of G if 
(a) K is compact and convex. 
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(b) K is the convex hull of K n 

(c) G is convexly linear on K, i.e. 

G(J2p^P^^ = if all Pj e K 

K is called complete, if all pure states which can appear in an optimal decomposition of any p E K 
are contained in K. 

Proposition 3.6 

With the assumptions of proposition 3.5 it holds: For g'~' (respectively g'^) and any cu E ^ there are 
complete convex leaves containing to. 

It makes sense to call the set of all complete convex leaves of G the convex foliation of G or, shortly, 
the G-foliation. 

Proof. The proposition will be proved for g'~'. The case of g'~^ is similar. There is an affine functions I 
such that 

l<g^, l{u;) = g''{uj) (35) 

We define a subset KofVlhy 

K = {p\l{p)=g''{p)} (36) 

□ 

Proposition 3.6.a: (36) is a complete g^-leat. 

Clearly, cu E K. Let us choose p E K. For an optimal pure decomposition, p = J^Pj'^^j^ obtains 

However, /(tTj) < g^{Trj) by assumption. By the equation above, all these inequalities must be 
equalities. Hence, the pure composers ttj of every optimal decomposition of any p E K are contained in 
K, i.e., K is complete. Now choose another p' E K and let p' = be an optimal decomposition. 

For < p < 1 we get, applying first our assumption and then convexity of g'~', 

l{pp + (1 - p)p') < g'^ipp + (1 - p)p') < pg'^ip) + (1 - p)g''{p') 

By the assumption the right hand side can be written 

plip) + il-p)l{p') = l{pp + il-p)p') 

and is equal to the left expression. Hence, equality must hold and K must be a convex set. We have 
seen already that every p E K can be represented by a convex combination of elements from K fl QP^'^^. 
Because g is continuous, the set of all tt satisfying g{Ti) = /(vr) is compact. Hence, K is compact. 
Indeed, it is convexly generated by a compact set. □ 

We repeat a further standard notation. An element p of a convex set K is called K -inner or "convexly 
inner" if for any u E K, and for small enough positive s, it follows (1 + s) p — sv E K. Geometrically, 
the line segment from u to p can be prolonged a bit without leaving K. There is also a topological 
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characterization: A i^-inner point is an inner point of K with respect to the affine space generated by K. 
(Example: The invertible density operators are the convexly inner points of Vl{1-L).) 

The intersection of complete convex leaves is either empty or it is a complete convex leaf. Hence 
there is a minimal complete convex leaf containing a given u E Vt, the co-leaf of g^. It is convexly 
generated by all those tt G iV^^"^ which can appear in an optimal decomposition of cu. 

The cj-leaf is the largest convex leaf containing cu as convexly inner point. 

Now let Ki be the wi-leaf and K2 that of U2. If oji is a convexly inner point of K2, then Ki C K2. In 
particular, Ki = K2 if and only if they contain a point which is commonly inner. Let us draw a corollary: 

If K2 is properly larger than Ki, the convex dimension of K2 must be strictly larger than that of Ki. 

A chain Ki C K2 C . . ., consisting of different complete (7^-leaves, cannot contain more than n + 1 
members. As above, n denotes the convex dimension of The maximal number of different leaves in 
any chain is called the depth of the ^(^-foliation. As said above, the ^(^-foliation consists of all complete 
(7^ -leaves. 

In the case = Vt^H), the depth if bounded by d?. If the roof is an affine function, see example 2.3, 
the faces of VL are exactly its leaves and the bound is reached. 

Let us look again to the setting above in more geometric terms. We shall see that the tj-leaves of 
and g'^ correspond uniquely to the faces of Vt[g]. 

The triple {^^[(7], 11} is di fiber bundle with bundle space VL[g], base space Vt, and projection 

n : {uj,\} ^uj (37) 
In this scheme a roof G becomes a cross section, say sg, by setting 

u^SG{uj) = {uj,G{uj)}en[g] (38) 

we get Yi{sG{oj)) = cu, which is necessary for a bundle structure. 
The boundary, of ri[g] is the union of three disjunct sets: 

d'Q[g] = {{co,g''{co)}\g''{co) = g''{cu)} (39) 
d-n[g] = {{cu,g''ico)}\g''ico)^g''icu)} (40) 
d^n[g] = {{uj,g''icu)}\g^ico)^g''icu)} (41) 

The cross section (38) with G = g^ maps Vl onto U d^^l[g] while the cross section with G = g^ 

maps the base space onto d'^Vl[g] U d^VL[g]. 

The fibres degenerate to a point at the boundary part (39) and can be identified with a subset of Vt. By 
proposition 3.2 that subset consists of the pure point and possibly of some faces at which g^ = g'^, and 
all roof extensions of g coincide and are affine. 

Let us consider a face K contained in the "lower" part d^Vt\g] of the boundary (40). K is convex by 
definition and compact because of the compactness of i7[(7]P"'^'^. Therefore, the projection K of K to Vl, 
n ^ = i^T, is convex and compact, see (37). We use the supposed face property: cu E K implies 

a) = '^^pjOJj =^ Uj E K (42) 

if all pj > 0. Writing this out in the manner u = {cu, g'~'(cu)}, and so on, we arrive at 

{u;,g''iu)} = 5^p,{a;„(7^(u;,)} = {a;, 5^^,(7^(0;,)} (43) 
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Now we can state 
Proposition 3.7 

There is a one-to-one correspondence between the faces of VL[g\ contained in yJd~VL[g\ and the 

complete convex leaves of g^. The cross section sg, with G = g^, maps complete convex leaves of g'~' 
onto faces of ^[g]. The bundle projection 11 returns them back to Vl. 
For the concave roof things are similar. 

3.3. Illustrating Examples 

Example 3.1: Minimum and maximum of g 

On let (jfmin be the minimum of g and g^^ax its maximum. The convex hull of the set 

{7ren^''''\gin)=g^,^} (44) 
is a complete convex leaf of g^. The convex hull of the set 

{7renP^'''\g{7r)=g^,.,^} (45) 

is a complete convex leaf of g'^. 

Let us again consider the diagonal map D{X) = diag(X) as in (16) and its von Neumann entropy 
S{D{.)), see (17). g{iT) = S{D{-k)) is the output entropy of the pure state tt. The Hilbert space 
dimension is denoted by d. As well known, the minimum output entropy is zero and the maximal 
one log d. 

Things become more refined by restricting the channel onto a face ofVl. As an example we take a short 
look at the {d — 1) -dimensional subspace T-Lq which is orthogonal to the vector \(p) = c?"^/^ Yl 
consists of vectors J2 ^jli) such that ^ aj = 0. "Hq supports some pure states satisfying D{7c) = d^^l 
and the maximal output entropy is log d again. 

There is a reasonable conjecture, saying that the minimal output entropy is independent of d and 
equal to log 2. There are d{d — l)/2 pure states tTj^, j < k. The matrix elements a„m of iijk are 1/2 for 
n = m = j and n = m = k. They are —1/2 for n = j,m = k and m = j,n = k, and all other entries 
are zeros. Hence it is evident that S{dia.g{TTjk)) = log 2 and, therefore, 

£;^(p)<iog2, p = ^(i-|^)(<^l) 

because we can represent p by a convex combination of the pure states Tijk. The conjecture asserts 
that the decomposition is an optimal one. The conjecture rests on the fact that for no other state than 
TTjk, supported by T-Lq, the output of the diagonal map is of rank two. Then one applies a theorem of 
Michelson and Jozsa, see appendix of [34], reducing in the case at hand the minimization of the entropy 
to that of minimizing the second elementary symmetric function or, equivalently, to the minimization of 
the concurrence. For d = 2 this is trivial, for c? = 3 it can be done, for d > 3 it is yet a conjecture. 

Under the assumption, the conjecture is true, the convex set generated by the vr^fc is a complete leaf. 
The minimal length of an optimal decomposition of p is equal to d{d — l)/2. 
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Example 3.2: Again the diagonal map 

We return to the example 2.4 to prove the flatness of the Entropy (17) of the diagonal map (16). Hence, 
since we know from von Neumann's work the concavity of — )■ S{D{ijj)), flatness of (17) proves the 
diagonal map a flat concave roof for all dimensions dofH. In other words, for 

g{7r) = 5(diag(7r)), vr G QiU) (46) 

it follows for any density operator u 

g'^iu) = 5(diag(a;)) (47) 

Proposition 3.8 

Choose uj e VliT-L). The set of all states p with diag(p) = diag(ci;) is a complete convex leaf of 
(47). g'^ is a flat concave roof. 

The set is compact, convex, and S{D{.)) is constant on it. Transversal to the sets the function 
S{D{.)) is strictly concave. This excludes that two density operators with different diagonals can belong 
to a convex leaf. Hence, the sets are convex leaves. 

4. Wootters' Method 

We are going to describe the fundamental idea in [13], see also [12], and its generalizations [14]. After 
a short introduction to anti-linearity, which is on the heart of the method, we present slightly simplified 
proofs for a class of convex and concave roofs. For reasons of uniqueness we sometimes write (., .) for 
the scalar product instead of Dirac's (. |.). 

The use of anti-linearity [35] goes back, at least in physics, to Wigner. He applied it to the time 
reversal symmetry [36] and he discovered the structure of anti-unitary operators, [37]. A highlight in 
the further development of this line of thinking is the proof of the CPT-theorem within Wightman's 
axiomatic quantum field theory. 

4.1. Anti-Linearity in Short 

We start with some elementary remarks. An anti-linear operator, say ^9, obeys the rule 

H ai\<pi) + a,\<P2) ) = ai^|0i) + a^^l^s) (48) 

An important fact follows immediately: Because of cd = 'dc* the eigenvalues of i!) form a set of 
circles. Indeed, if \x) is an eigenvector of 13 with eigenvalue a, then €\x), \e\ = 1, is an eigenvector 
with eigenvalue e*a. Consequently, most of the unitary invariants of linear operators are undefined for 
anti-linear ones. The trace, for example, does not exist for anti-linear operators. 
The Hermitian adjoint i!)'^ of an anti-linear operator i!) is defined by 

(01,^^02) = (02, (49) 

There is to set a caution mark: Do not apply an anti-linear operator to a bra in the usual Dirac manner! 
By (49) one may get absurd results. 
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A useful class of anti-linear operators are the Hermitian ones. By (49) any matrix representation must 
result in a complex symmetric matrix. About symmetric matrices see [38]. It follows that the Hermitian 
anti-linear operators constitute a complex linear space of dimension d{d + l)/2 if dim'H = d. 

An anti-unitary, V , is an anti-linear operator which is unitary, i.e., satisfies = V~^. A conjugation, 
9, is an anti-unitary which is Hermitian. It implies 6^ = 1. In accordance with what has been said about 
eigenvalues, one can find an orthogonal basis . . . such that 6|0j) = with arbitrarily chosen 

unimodular numbers ej. 

A conjugation 6 distinguishes a real Hilbert subspace "He of H consisting of all 6-real vectors, 

There is a polar decomposition, = for any anti-linear operator \-t}\ denotes the positive root 
(^t^)i/2 y anti-unitary operator. The proof is similar to the linear case [39]. 

Now we turn to the case of an anti-linear Hermitian operator 'd = 'd'^ . It commutes with the positive 
(linear!) operator 'd'^ and, therefore, with I'd] = (^9^)^/^. With a non-singular i!) we perform = 
the square of which is 1. As it is Hermitian too, it is a conjugation. We conclude 

^ = e|?9| = |?9|e, 6 = 6^ = 6-^ (50) 

By continuity, or by a more detailed analysis, (50) can be verified for all anti-linear Hermitian i!). 

4.2. Building Roofs with an Anti-Linear Hermitian 'd 

Let be Hermitian and anti-linear on a Hilbert space V. of dimension d. This setting provides 
a function 

g{7r) = \{^P,^^P)\, 7r = |^)(V^| (51) 
on the pure states of fi('H). Now g'^ and g'~' are well defined and we shall prove: 

Proposition 4.1 

Let g be as in (51). If {Ai > A2 > . . .} denote the eigenvalues of \^/cu'^^/Lo\, then 

g'^iu) = max{0, Ai-^A,}, (52) 

g'^iu) = (53) 

g'^ and (7^ are flat roofs. 

At first we simplify the assertion by starting with 

t9^:=v^?9v^, \^J = i^uMV^f^ (54) 

Up to normalization every pure decomposition of u can be gained from a decomposition of the unit 
operator 1, 

U = ^ a/cJvTj a/w, 1 = ^ TTj (55) 

(51) is 1-homogeneous on the positive rank one operators. Comparing (55) with (25) and (26), it can be 
seen that 

= inf5^|(V'„^.V',)l = sup^|(V'„?9.^,)| (56) 
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where we have to run through all rank one decompositions of 1, 

Now "d^ is Hermitian and anti-linear. Therefore there is to any chosen set of phase factors ei, e2, . . . a 
basis ipi, (y92, • • • satisfying 

l^<^l Wj) = Wj) = (58) 

and the conjugation in the polar decomposition multiplies the j-th basis vector by tj. 

For the next step we assume the existence of a real d x d Hadamard matrix. Then we can choose a 
basis fulfilling 



\Xi 

\/d ' . 

3 

,2 



because of the orthogonahty and a^j = 1 we get for all k 

d{Xk,'&i^Xk) = ^akiakj{^i,'&i^^j) = ^CjAj (60) 



Therefore, by (51), we get 

9''{uj)<\J2e^X,\<g''{uj) (61) 
By varying the unimodular numbers ej, which could be chosen arbitrarily, one arrives at 

g'^iuj) < max{0, Xi-J^ > J2 ^^2) 

Assuming equality in (62), we see from (60) that a; is a flat point of g'~^ and of g^. If g'~'{co) > 0, then 
we choose ci = 1 and ej = —1 for j > 1 to obtain from (58) an optimal basis {|xfe}- In the concave case 
we set ej = 1 for all j. 

Now we are going to prove equality in (62), starting with g'^. We clearly get an estimation from below 
in (56) by 

inf I , = {i^k,'&uj'4'k)\ (63) 

Sandwiching with the eigenbasis (58) it yields 

X^^j = ^{-^k^^Lo^j) {^j,^k) = ^ejXj{tjjk,(pj) {^j.tpk) 

jk jk 

and this shows, summing first over k, 

k k 

Its minimum is attained by the largest of the two numbers and Ai — A2 — . . . and the first equation in 
(54) is true. 

Now we prove equality in (62) for g'^. the second equation of (56). It is 

\{^k,^u.i^k)\ = \{^k,Kmi Wk) = Q\^k) 
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We apply Cauchy's inequality. The result is 

and, because O is an involution, hence anti-unitary, we arrive at 

\{il^k.^M\ < {ipkMu.\ilJk) (64) 

Summing up we get the trace of 1^9^^! which upper bounds g'^{uj). 

Up to now the proof of proposition 4.1 rests on particular bases {|xfc)}- They exist if there is a real 
d X d Hadamard matrix, d = dimH. To overcome the restriction we go to a larger Hilbert space, Ti^T-Lo, 
for which the proposition has been proved. Then we restrict to the face of density operators supported 
by the original Hilbert space T-L. Because a (flat) roof remains a (flat) roof if restricted to a face, the proof 
will become complete. 

To do so, we choose d' = dim?/ © "Ho sufficiently large and extend to 'd' by requiring 'd'l'ipo) = 
for all iV^o) £ ^0- We then choose any conjugation ©o on T-Lq and use 6' = O © Qq. Now, if there is a 
real d' x d' Hadamard matrix, we are done. 

Hadamard matrices exist in dimensions rf' = 2™. This suffices for the proof. 

Proposition 4.2 

Let g{n) = \{'ip,'dip)\. Then and g'^ allow for flat optimal decompositions of length d' where 
d < d' = 2^"-. More generally, if d < d' and there is a real d' x d' Hadamard matrix, then there are flat 
optimal decompositions of length d'. 

4.3. Cases of Application 

Indeed, the question is now: How to find a suitable anti-linear Hermitian operator i!) to calculate 
concurrence, tangle, and entanglement entropy (as particular cases of entanglement of formation) 
in 2 X n systems. Clearly, this can be fully successful for flat roofs only. 

Let T be a trace preserving positive map from the states 

u en'^ ■.= n{n), dimn = d 

into the 1-qubit state space Vl{1-L2)- Because the trace of the output, Tr T{uj) is one, T{uj) is characterized, 
up to a unitary transformation, by one variable. It is common to use 4 det Tiu) or its square root to be 
this variable. Let us abbreviate the convex roofs on Vf^, playing a role below. By 

Ct = 2(Vdetr)^, (detT)(w) = detT(u;) (65) 
Tt = 4(detT)^ (66) 
Et = (5t)^ Sriou) = S{T{u)) (67) 

(65) is the concurrence, (66) the 1 -tangle, and (67) the entanglement of T. (67) is the entanglement of 
formation if T is a partial trace of a bipartite quantum system. By its very definition we need only the 
values of det T and of St for pure input states. 
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There are some general relation between these three roofs. The first is typical also for more general 
settings: Ct is a positive convex function and so does its square. For pure input states the tangle and 
the squared concurrence coincide. Hence, because r-p is maximal within all convex extensions, it is not 
less than C|.. On the other hand, if u turns out to be a flat point of Ct, than this remains true for its 
square. Thus, 

Proposition 4.3 

For stochastic maps with 1-qubit outputs it holds 

tt{uj) > CricoY (68) 

and equality takes place for flat roof points of Ct- 

Let us now switch to St- With use the abbreviations 7]{x) = —xlogx and 

ax) = J^-^]+J^-p-], l = x' + y' (69) 



^ is defined and continuous for — 1 < x < 1 and it is strictly convex. Therefore, ^ is the sup of a family 
of functions ax + b. Inserting a convex function C defined on any convex set with values — 1 < C < 1, 
we get ^{C) by a sup of convex functions aC + b. Therefore, ^{C) is a convex function on the domain 
of definition of C. 

We apply this fact to the concurrence Ct yielding: 

w^^tM :=e(CTM) (70) 
is a convex function on ^l'^. Ct for pure states tt we get 2A/detT(7r). We insert in (70), 



w) = vi — - — ^ ) + vi — - — ^ ) 

and one identifies the arguments in r] as the two eigenvalues of T(7r). 

^tW = 5(T(7r)), TTGfi'^'P"^^ (71) 

proves ^T to be a convex extension of the pure states output entropies. Reasoning as for proposition 4.3 
results in 

Proposition 4.4 

For stochastic maps with 1-qubit outputs it holds 

^tM >e(CTM) (72) 

Equality takes place if w is a flat roof point. 

It should be underlined hat there are more and different estimations for concurrence and entanglement 
of formation in higher dimensions, see [40-44]. In [56] there is an application to states with only two 
different, but arbitrarily degenerated eigenvalues. 
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4.4. How to Find 

There is a general recipe to get the wanted anti-linear operator for channels T mapping the states of a 
quantum system Ti^ to 1-qubit states. Assume a Kraus representation 

T(X) = ^ Aj : nd^n2 (73) 

There is an additional condition to be fulfilled: T must be Kraus representable with not more than two 
Kraus operators. But at first we remain within the more general (73). 

The key to the following is the existence of the "time reversal" or "spin flip" anti-unitary operator Of 
on 7^2, 

ef(co|0) + ci|l)) = ct|0)-c*|l) (74) 

Apart from the obvious 

$1 = 6f ^ = — 6{ 

the relation 

OfY^OiY = -{detY)l2 (75) 

is valid for all operators Y E B{'H2)- Up to a multiplicative constant, only the spin flip commutes with 
all U E SU(2), It is really a very special anti-unitary operator. 

The task is in inserting Y = T{X) into (75) and to get something similar for det T. This goes through 
particulary nice if X is of rank one, X = \-ip2) (V'll- Calculations show 

detT(|V^2)(V^i|) = Yl {^2,^jk^2r, \^^) e (76) 

j<k 

The anti-linear Hermitian operators 'djk are defined by 

^,k = l{A*OiAk-AiefA,) (77) 

using the Kraus operators Aj from any Kraus representation (73) of T. The operators -^jk are Hermitian 
and anti-linear. 

In the lucky case of channels (73) with just two Kraus operators, Ai and A2, we get only one operator 
1? by (75) and, hence, 

VdetT(|V)(V'|) = |(V^,^V^)| (78) 

and, by proposition 4.1, we are done. 

Wether and how one can replace the operation X — > 'dX'd by an anti-linear stochastic map to obtain a 
more general roof construction, is unknown. An (implicit) attempt is in [46] by taking X — )• (Tr X)l—X 
as an higher dimensional substitute for the flip operation. 

4.5. Applications 

The partial trace of a 2-quibt system can be represented by two Kraus operators: Looking at the 
operators over ® I-L2 as block matrices, the partial trace over the second part is the map 

X=[^l -^Xoo + Xn = F (79) 
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One immediately sees a possible choice for the Kraus operators, 

= -^{l2, O2}, A2 = -^{02, I2} (80) 
Now we can compute d according to (77). We get, eventually up to a sign, Wootters' conjugation, 

^f' 



4^9 = ^„=(^_^^ ^j=e,®e, (81) 

Remark: The concurrence of a 2-qubit system is a flat convex roof. Generally, complete convex leaves 
consists of a set of flat ones. 

Let us choose two operators, Ai = A, A2 = B, from B{'H2) so that (73) becomes a 1-qubit channel. 
Then there are only two eigenvalues, Ai, A2, of \^/LO'^^/uJ\ to respect and (52) simplifies to 



g'^ico) = IA2-A1I (82) 
One has to solve a quadratic equation to get the general expression 

^ Cricof = Tr (coM) - 2(det X) (det ^^^^^ (83) 

There are standard forms for 1-qubit channels, [20,47^9]. For channels with two Kraus operators 
one can assume 

A = i , 5 = , (84) 

up to unitary equivalence. With this choice i!) acts as 

i9(co|0) + ci|l)) = (6ioaooCo)lO) - (6oiaiiCi)*|l) (85) 

After inserting in the relevant expressions one observes that the concurrence is the restriction of a 
semi-norm to the state space. Indeed, one obtains, [50,51], 

Ct{X) = I Ifeioaoolxoo + l&oicnikii + zxw - z*xoi \ (86) 

z is one of the roots of 

z^ = alQauboiblo (87) 

5. A Subtraction Procedure 

We start with a simple example not covered by Wootters' method and showing the principle. The idea 
is to subtract from det T{X) a suitable multiple of det{X) to get the squared concurrence or the tangle 
of a stochastic 1-qubit map T. We choose for T the map 



3^00 2:01 \ ^ / xoo + (1 - 7)2^11 

xio xuj V 7x11 



(88) 



and consider 

7x00X11 + 7(1 - 7)xn - w(xooa;ii - xoiXio) (89) 
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With = 7, we arrive at the squared concurrence of T, 

^ Cr(X)2 = 7(1 - 7)x?i + 7x01x10 (90) 

At first (90) is a positive semi-definite quadratic form implying that its square root is convex. Secondly, 
by its very construction it coincides for pure states vr with 4 det Tiji). Finally, on the state space, it is a 
roof. To indicate a general way to proof the roof property we polarize (90) and get 

(X, Y)t := 7(1 - 7)xiiyii + ]p{xQiyio + xio^oi) (91) 

which is a positive semi-definite bilinear form in the space of Hermitian matrices. For the pure state 
TTo = |0)(0| (90) becomes zero. Applying the Schwarz inequality we get (X, ttq) = for all X. Hence, 
with any pure state vr, 

ujs = {I - s)7ro + STT ^> [us, uJs) = s'^{tt, tt) (92) 

Taking the root we see that s — > Us, < s < 1, is a convex leaf of Ct for every vr. Now the assertion is 
proved and, the more, Ct is not a flat roof. 

One may ask, whether one can diminish w a bit without destroying convexity on Vl. Let us use in (89) 
the new value w' = 7^. Then (89) becomes 

7x00X11 + 7(1 - 7)xii - 7^(xooXii - xoixio) 

Different to the former case we restrict ourselves to ^2(7/2) and respect the condition xqo + Xu = 1. 
After some manipulations we arrive at a convex roof which coincides with det T(7r) for pure states. Up 
to a factor it must be the tangle of T on fi. 

|rT(X) = 7(1 -7)xii + 7^X01X10 (93) 

The tangle is affine on the set of density operators with xoi = constant. As the square of the concurrence 
is equal to the tangle for pure states, we have the inequality tt{cu) > Ct(c<;)^ for mixed states. 

That the ansatz (89) is working generally for the concurrence has been shown first in [16,17,52], by 
means of the "S-lemma of Yakobovich" and, using the explicit expression for general stochastic 1-qubit 
maps of [20], in [19]. The case of the tangle can be read off from [22]. 

5.1. Concurrence of Stochastic 1-Qubit Maps 

Let T be a stochastic, i.e., a positive trace preserving map. We prove 

Proposition 5.1 

There is a real number < w < 1 such that for all p G ^(^^2) 

1 Cripf = det T(p) - w det p (94) 

At first we show the w-bounds. With w < 0, (94) becomes the sum of two concave functions on 
and Ct could not be convex. To prove 1 > w we insert po = (1/2)1 and get det T(po) > w/A. However, 
det T(po) < 1/4 is required by stochasticity. w > 1 would be a contradiction. 
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Next, we consider the expression (94) on the Bloch space of all Hermitian operators of trace one, 

X = ^(1 + Xi(Ti + .X2Cr2 + XsCTs) (95) 

The determinant of X is a quadratic function 

det X = -^{1 — x\ — x\ — x\) 

in the Bloch coordinates. The same is with det T{X). In the terminology of the S-lemma, the quadratic 
function detT(X) is called co-positive with detX because from detX > it follows detT(X) > 
under the constraint Tr X = 1. 

The S-lemma, [18], states the following: Let qi and q-z be two quadratic functions on WC", not 
necessarily homogeneous. If qi{x) > implies q2ix) > then q2 is co-positive with qi. 

S-lemma 

Let qi be strictly positive at least at one point. Then there exists a real number w such that 

q2{x) - wqi{x) > for all x G IR" (96) 

if and only if q2 is co-positive with qi. 

Now it is proved: The right expression of (94) can be made non-negative with suitable numbers w for 
all Hermitian trace one operators. Then, substituting Xj — )■ xj/xq and multiplying by x^, we see that 

det Tr(X) - w det X > (97) 

becomes a positive semi-definite and homogeneous polynomial in the variables Xq, Xi,X2, x^ for some 
w. Polarizing (97) as in (91) yields a positive semi-definite symmetric form (X, Y)^ satisfying 

(X,X)^ = detT(X) - wdetX (98) 

We now may assume (97) for all values w which are bounded hy 1 > wi > w > W2 > 0- We can further 
assume that (X, Y)^ is degenerated for w = wi and for w = W2. But wether degenerate or not, (97) 
implies Cauchy's inequality 

\{X,YW<{X,XUY,Y)^ (99) 

In particular, if F = is in the null-space, {u, v)^, = 0, then det T{v) = w det u and det T(z/) > w' det u 
for all allowed w', i.e., 

{w -w')detu >0, Wi>w'>W2 (100) 

Furthermore, (z/, X)^ = for all Hermitian X. 

We like to show: lfw = W2 every density operator p is a roof point. 

We have to distinguish several cases. The first one is Tr z/ = 0. Then, with p G fi, we define 

Ps=P + su (101) 

On this line {ps, ps)w is independent of s. Therefore, because of (94), its intersection with 1] is a flat 
convex leaf of Ct. 
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Being Hermitian and not identical zero, Tr z/ = results in det z/ < 0. Because of (100) one concludes 

For the next cases we suppose Ti v = 1 and start with 

ps = {l- s)iy + sp (102) 

Then {ps, Ps)w = s^- If the sign of s does not change while ps is inside the Bloch ball, we can take the 
root and get a convex leaf. This takes place if u is not an inner point of Vt, i.e., not a properly mixed state. 

The condition is certainly satisfied if det z/ = and z/ is a pure state. The condition implies the 
w-independence of (z/, z/)^ which becomes equal to det T(z/). As there is only one convex roof Ct it is 
wi = W2 necessarily [23]. 

u is outside if det u ^ 0. We conclude w = W2 from (100) and we get convex leaves by intersecting 
the lines (102) with the Bloch ball. 

In addition one observes: det z/ 7^ is necessary for wi > W2. Indeed, by (100) we see (z/, z/)^ = 
can be satisfied if either w = Wi and det u > orhy w = W2 and det v < 0. 

Now, all relevant cases are discussed and proposition 5.1 is proved. 

In the course of the proof one obtains two more general insights: 

Proposition 5.2 

The concurrence of stochastic 1-qubit maps is the restriction of a Hilbert semi-norm to the state space. 
Every state allows for an optimal decomposition of length two. 

5.2. Axial Symmetric Maps, Concurrence 

We shall list the subtraction parameter w for the class of axial symmetric stochastic maps. A standard 
form for them reads 

_ / axQQ + (1 - 7)3^11 Pxoi 
V /3xio 73^11 + (1 - 

with real non-negative parameters a, /3, 7. The trace preserving is obvious. Positivity requires < a < 
1, < 7 < 1, and 

< /^Lx := 1 + 2a7 - a - 7 + 2^ a{l - a)^{l - ^) (104) 
T is a channel, hence completely positive, if < 0^7. To express w one needs the "critical" (3c 

:= 1 + 2a7 - a - 7 - 2^a{l - a)-f{l - 7) (105) 

Then, see [52], 

w = max{l3^,l3'^} (106) 

At the bifurcation point (3 = (3c the T-concurrence is affine on Vt. For (3 > (3c the roof is flat. 
Otherwise it looks similar to the particular one (88). See [52] for more details. 



(103) 
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5.3. Axial Symmetric Maps, Tangle 

The tangle for stochastic 1-qubit maps can be found in [16,17]. For 1-qubit channels it is already 
in [22]. These tangles always allow for optimal decompositions of length two. 

The axial symmetric maps (103) can be treated explicitly. The following is due to [53]. 
Case A: If > |q; + 7 — 1| one has to use w = (3'^. It results in 

^TriX) = 1 - /32 _ (« _ _ 2(« - 7)(a + 7 - l)x3 + [/^^ - (« - 7 - l)']^^ (107) 

Case B: Here |/3| = |a + 7 — l|,a bifurcation point in the parameter space, w = (3'^ results in 

^MX) = 1 - - (a - 7)2 _ 2(a - 7)(a + 7 - 1)0:3 (108) 

and the tangle becomes affine on the Bloch ball. 

Case C: If < |a + 7 — 1| then to = (a + 7 — 1)^ and we obtain 

^tt{X) = 1 - (a + 7 - 1)2 - (a - /3)2 - 2(a - 7) (a + 7 - l)x3 + [(a + 7 - 1)' - {4 + xl) (109) 
6. Symmetries 

The use of symmetries is almost obligatory in the treatment of roofs. We present only a small, 
hopefully helpful, part of it, mainly abstracted from [10,26,27,54,55]. See also [45]. 

Let fi('H) be the space of states supported by the Hilbert space T-L of dimension d, and g a real 
continuous function on Q'^^^'^. 

A symmetry of i7 is a transformation 

Lo ^ := VooV^^ (110) 

is a unitary or an anti-unitary operator inducing the symmetry (110). 
We need the group T of all V such that 

vrGfiP-- ^ g{n'')=gin) (111) 

r is the invariance group of g. A quite obvious statement reads 

Proposition 6.1: 

If r is the invariance group of g, then 

^7^(a;^)=^7^(a;), g^iu"") = g'^iuj) (112) 

for all V e r and all u eVt. 

Let iiT be a convex leaf of g^. Obviously, 

K^ = {uj'^\ueK} (113) 

is a convex leaf of (7^ again. Hence, K — )■ permutes the convex leaves. 
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An interesting subgroup of F is the stabilizer group 

Vk = {V eV\K^ = K} (114) 

This is a compact group with an invariant Haar measure, say cLkV . We can perform the invariant 
integration ("twirling") over Vk, 

u^u^ = I u^dxV (115) 



There is only one Fj^-invariant element in the convex hull of all tu^ . It is co^ . 
The map u — t- uj^ contracts K onto the set 



^stable = {^^K\U = J U^dxY} 



(116) 



This set, being convex and compact, is the convex hull of its extremal invariant states. Extremal invariant 
states can be represented by vr^ with pure n E K. Invariant states which are not extremal, cannot be 
represented in such a way. 

Proposition 6.2: 

Let K he a convex leaf of and F^ its stabilizer group. Then K fl f2P"'^'^ consists of F/^-orbits. 

Every extremal Fj^-invariant state of K is of the form vr^, vr pure. 

Every Fi^-invariant states of is a convex combination of extremal FA -invariant states. 

6. 1. Entanglement of the Diagonal Channel 

In example 3.2 we considered the concave roof of (46), 

g{Ti) = 5(diag(7r)), vr G ^iU) 

Now we look at the convex roof g^, i.e., at the entanglement Ed of Diu) = diag{uj). 

With the exception of d = 2 one does not know the structure of Ej:,. But there are some insights on 
highly symmetric, "isotropic" quantum states. 

The channel D will be described by the help of a basis j = 1, . . . , ci. To it we associate the vector 



y/d 

The invariance group, F, of (1 17) consists of the permutations of the chosen basis, eventually followed 
by the conjugation defined by 0|j) = |j) for all j. The density operator tu commutes with F if its 
matrix elements satisfy 

{jHj) = 1, (Jk|fc) = ^ (118) 
for all j and all /c 7^ j. x is a real number in the range 

--^<a;<l (119) 
d — 1 
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The restriction is due to the positivity of u. One often uses the fidelity parameter, F, 

< F := {^Icol^) = < 1 (120) 

and we denote the corresponding F-invariant density operator hy up- 

We choose an allowed value F and write K^^ or K(F) for the complete convex leaf ofcop with respect 
to Ed- Now we state the following: 

pe ^ ep = pB (121) 

6 is the conjugation about the basis 

Proof: It suffices to prove the assertion for pure states. We assume that the pure states vr and vr' = Ovr© 
are both in K^^. The diagonal parts of them and of p = (l/2)(7r + n') are the same. It follows Snip) = 
Soi^^) = 80(7^') and, because we are inside a convex leaf, we get also g'~'{p) = Sd{p)- Hence, by 
proposition 3.2, Sd = Ed on the whole face containing vr and vr'. This is a contradiction if vr 7^ vr'. 
Indeed, we proved something more: 

Proposition 6.3: 

Any two different pure states contained in a convex leaf of E^ must have different diagonal parts. 
The F-invariant density operators are ordered by the fidelity parameters as indicated by (120). The 
more, F ^ cupis convexly linear in F. This fact allows to apply proposition 6.2. 

Proposition 6.4: 

Let K(F) denote the compact convex leaf of Ed belonging to the F-invariant density operator tup- 
There is a subset R of the unit interval as follows: 

a) Either we have F E R. Then K{F) consists of flat roof points only. There is no other maximal 
symmetrical state in K{F) than cop. 

b) Or there are F", F+ e such that K{F) is the convex hull of K{F-) f] K{F+). 
In the case c? = 3 and if 

1 8 

_ < ^ < p** p** = _ 

2 - - ' 9 

or F = case a) is true. There is an optimal vector of the form a\l) + b\2) + b\3) for any of these 
F-values. K^^ = K{F) contains not more than three pure states. They become permuted by the action 
of the invariance group F. 

For more details, also in higher dimensions, see [55]. 

6.2. An Embedding 

There are numerous relations between different channels. Some of them provide insight in roof 
structures. Our main interest is again in the diagonal channels. 

Let H and T-L' be Hilbert spaces of dimensions d and d' > d. There are embeddings of T-L into T-L' 
which relate the entanglement Ed and Ed' of the corresponding diagonal maps, D and D'. 
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For their description we first choose d integers, mi, ... , m^, such that d! = mi + . . . + m^ and 
enumerate the basis vectors of T-L' as 

\jk), j = l,...,ci, k = l,...,mj (122) 

We further choose numbers 

yj,k, j = l,...,rf, k = l,...,mj (123) 

satisfying 

J2\yjk\' = ^, 3 = l,---,d (124) 

fc=i 

These data provide a unitary embedding 

d 

\3)^V\3) = Y,yjk\3.k) (125) 

From 

diag(X) = {xii, . . . , Xdd] (126) 

we get 

diag(FX\/"^) = {lyiipxii, . . . , lyimj^xn, \y2i\^X22, ■ ■ ■ \y2mA^X22, ■ ■ •} (127) 
One obtains for the entropy of the diagonal channel 

mi m2 

J2vi\y2k\'x22) + . . . (128) 

fc=l A;=l 

The functional equation r]{xy) = yr]{x) + xr]{y) allows to rewrite the first sum in (128) into the form 

\yik\'^r]{xn) + xn^ridyikl'^) 

and we get finally 

So'{VuV^) = SD{oo) + Y.{Mj)Y,vi\yjk\') (129) 

j=l k=l 

A convex roof remains a convex roof if we add a function linear in ui. Plugging (129) into (25) 
directly provides: 

Proposition 6.5: 

If y, is embedded unitarily in T-L' according to (125), then the entanglements of the diagonal channels 
are related by 

ED'iVioV^) = Ed{uj) + /(w) (130) 



and the linear function is given by 



1{^) = Y.^3W\3)Y.'^M) (131) 

i=l k=l 
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Optimal decompositions are mapped onto optimal decomposition and convex leaves onto convex leaves. 
A particular simple example is the embedding 

\/|0) = |l), \/|l) = -L(|2) + |3)) (132) 

of I-L2 into 7/3. Pairs of pure state vectors, yielding flat optimal decompositions for the entanglement of 
D2, are 

ao|0) + ai|l), a\\Q)+al\l) (133) 



A particular case is Oq = i/l/3, ai = a/2/3. Applying the map (132) results in the optimal pair 



^(|1) + |2) + |3)), V^|l) + V^(|2) + |3)) (134) 
CaU TTo and tti the pure states determined by the vectors (134). Then 

diag(7ro) = {J, J, diag(7ri) = ^, ^} (135) 
3 3 3 3 o o 

and the fidelity parameters (120) are 1 and F** = 8/9. Returning to the remarks after proposition 6.4, 
we see a reason why this value should be a bifurcation point for the behavior of maximal symmetric 
states and their convex leaves. 

6.3. A Further Embedding 

There are strong relations between the entanglement of the diagonal channels and the entanglement 
of formation, governed by embedding procedures. A nice and quite simple one is 

V\3) = = \j) ® |j) (136) 

is a unitary map of "H onto the subspace H' of "H" ® V.^ with basis The embedding is of interest 

because of 

D{X) := diag(X) = TuiVXV^), X G B{n) (137) 

The relation implies 

E{{VioV^)) = Ed{uj), uj G n{n) (138) 

and relates the entanglement Ed of a diagonal map to the entanglement of formation E for bipartite states 
supported by T-L'. This is true for all dimensions d = dimH. The vector (117) becomes a completely 
entangled one, say |e), if transformed by (136). 

The crucial point is now that invariance group of |e) = V\tp) is much larger than VFV^ the 
invariance group of lip)- Tg contains the local unitary operators U ® U and the swap operation. The 
involution VQV^ = 6 (g) 6 is defined originally only on 1-L' . The canonical extension to an involution 
0e of all "H" ® 1-0' can be gained by 

Qe\3k) = \kj) (139) 
and the requirement of anti-linearity. 6e satisfies 

Q,{X®l)\e) = {X^ ®l)\e) (140) 
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for all X G Bi^H") [57]. 

Now one tries to enlarge convex leaves in 1-L' by transforming them with operators from a suitable 
larger group F' C Fg. The involution ©g is a symmetry of the entanglement of formation. If a bipartite 
state p commutes with ©e, all elements of the convex leaf of p must commute with ©e. 

One can understand quite well why Ed for maximal symmetric states are so similar to the 
entanglement of formation for isotropic states. A detailed discussion is not in the frame of the present 
paper. However, an essential point in the considerations above is in the relation between the entanglement 
of diagonal channels and the entanglement of formation. This may be of use in future research: It seems 
easier to imagine the structure of the diagonal channel as that of the partial trace. Nevertheless, the 
degree of difficulty is about the same. 

7. Summary and Outlook 

Given values g{'n) for pure states tt, the direct way of solving the convex roof problem is the search 
for optimal decompositions. The most prominent and successful examples are the concurrence and the 
entanglement of formation for 2-qubit bipartite quantum systems. The method goes back to Wootters 
and is described in section 4. With it one gets analytical expressions and flat optimal decompositions. 
The flatness of the convex (and concave) roofs inherited from Wootters' method is rendering its use in 
higher dimensions. 

A quite different way is to look for a maximal convex extension G of (7 : If for any other extension G' 
of g we find > G{ijj) for a state u eVL, then G' cannot be convex. 

A further, and more efficient reformulation of the convex roof problem asks for roof points, (see 
definition 2.1a), of a convex extension G of g. At a roof point w of a convex extension G one gets 
G{ijj) = g^{uj) and the problem is solved for the particular state oj. Similarly, if G would be a concave 
extension of g, then G{ijj) = g'^{uj) for a roof point of G. 

This way of proving is used in the chapter on the "subtraction procedure". One of its merits is the 
control on the concurrence and on the 1 -tangle of any rank two density operator of a 2 x m bipartite 
quantum system. The same is with the slightly more general class of stochastic, (just positive and trace 
preserving), 1-qubit channels. 

Therefore, there is some hope to get the concurrence (and the 1 -tangle) for all states of any 2 x m 
bipartite quantum system explicitly. 

But even if this becomes true, it does not provide us with the entanglement of formation of a 2 x 3 
system: The concurrence ceases to be flat. However, by proposition 4.1 one can obtain reasonable 
lower bounds. 

Wootters' and the subtraction method seem to be quite different in spirit. Uniting the strength of both, 
would be very useful. Also one should look at the subtraction method in more general terms. One can 
get at least lower bounds on the concurrence for higher dimensional system as shown by two examples 
in [52]. A more systematic study of the issue seems prospective. 

The use of symmetries is well know and efficient in general. In convex roof construction the 
symmetries of (7 as of a function on the pure states is what counts. If a state to is invariant with respect to a 
symmetry group F, its convex leaf, (see definitions 3.2), is the convex hull of a set of F-orbits consisting 
of pure states. The shapes of the leaves can be quite different. However, one would suppose a smooth 
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change of the leaves with the exception of some bifurcation points (or lines ...) at which the dimension 
of the leaves is changing. Some help comes from embedding a lower dimensional problem into a higher 
dimensional one. This is shown for the entanglement of the diagonal channel: Et can be computed 
in any dimension on 2-dimensional subspaces which contain at least one pure diagonal state. On these 
subspaces Et is the sum of a flat convex roof and a linear function. The study of more examples is 
certainly desirable. 

At this point we can return to the concurrence of the stochastic 1-qubit maps. For them Ct is the 
restriction of a Hilbertian semi-norm to the state space. The proof of proposition 5. 1 shows a one to 
one correspondence between the structure of the foliation and the null-space of the semi-norm. Two 
stochastic 1-qubit maps come with the same pattern of their convex leaves if the null-spaces of their 
semi-norms are identical. Indeed, it is the first class of channels and roofs with a complete classification 
of their convex leaves and their foUation. 
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